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We study a model where two scalar fields, that are subdominant during inflation, decay into 
radiation some time after inflation has ended but before primordial nucleosynthesis. Perturba- 
tions of these two curvaton fields can be responsible for the primordial curvature perturbation. 
We write down the full non-linear equations that relate the primordial perturbation to the 
curvaton perturbations on large scales, and solve them in a sudden-decay approximation. We 
calculate the power spectrum of the primordial perturbation, and finally go to second order to 
find the non-linearity parameter, /nl- Not surprisingly, we find large positive values of /nl if 
the energy densities of the curvatons are sub-dominant when they decay, as in the single curva- 
ton case. But we also find a novel effect, which can be present only in multi-curvaton models: 
/nl becomes large even if the curvatons dominate the total energy density in the case when 
the inhomogeneous radiation produced by the first curvaton decay is diluted by the decay of a 
second nearly homogeneous curvaton. The minimum value min(/NL) = — 5/4 which we find is 
the same as in the single-curvaton case. Using (non-)Gaussianity observations, Planck can be 
able to distinguish between single-field inflation and curvaton model. Hence it is important 
to derive theoretical predictions for curvaton model. From particle physics point of view it is 
more natural to assume multiple scalar fields (rather than just one "curvaton" in addition to 
inflaton). Our work updates the theoretical predictions of curvaton model to this case. 



1 Introduction to curvaton model 

Theories beyond the standard model often contain a large number of scalar fields in addition 
to the standard-model fields. In the very early universe it is natural to expect the initial values 
of these fields to be displaced from the minimum of their potential. If they are displaced by 
more than the Planck scale, then they can drive a period of inflation. Otherwise they will 
oscillate about the minimum of their potential once the Hubble rate, H, drops below their 
effective mass. An oscillating massive field has the equation of state (averaged over several 
oscillations) of a pressureless fluid. Thus the energy density of a weakly interacting massive 
field tends to grow relative to radiation in the early universe. Such fields must therefore decay 
before the primordial nucleosynthesis era to avoid spoiling the standard, successful hot big bang 
model. If the energy density of a late-decaying scalar is non-negligible when it decays, then any 
inhomogeneity in its energy density will be transfered to the primordial radiation^ Therefore, 
the observed perturbations in the Universe can solely result from these inhomogeneities and the 
perturbations in the inflaton decay products (i.e., pre-existing radiation) can be negligible. This 
is the curvaton scenario for the origin of structured We focus on this pure curvaton model@ 

2 Definitions: primordial perturbation, power spectrum, and /nl 

The primordial density perturbation can be described in terms of the non-linearly perturbed 
expansion on uniform-density hypersurfaces^ 1 



a Before their decay, the curvatons represent pure isocurvature degrees of freedom which are ruled out by data^ 
Nevertheless, the pure curvaton model generates adiabatic primordial density perturbations, if all the species are 
in thermal equilibrium and the baryon asymmetry is generated after the curvatons decay. However, if inflaton 
perturbation was not fully negligible and the curvaton(s) decayed into radiation and cold dark mater (CDM), 
then a correlated mixture of adiabatic and CDM isocurvature primordial perturbations would follow. These 
possibilities have been studied extensively, e.g., in Refs.^1 It should be noted that an admixture of subdominant 
isocurvature (up to 5-10%) and predominant adiabatic primordial perturbation can fit well^the current data. 




(1) 



where N = J Hdt is the integrated local expansion, p the local density and p the local pressure, 
and p is the homogeneous density in the background model. We will expand the curvature per- 
turbation at each order (n) as C(^ x ) = J2^=i ^C(n)(^ x ) > where we assume that the first-order 
perturbation, £n)> is Gaussian as it is proportional to the initial Gaussian field perturbations. 
Higher-order terms describe the non-Gaussianity of the full non-linear £. 

Working in terms of the Fourier transform of £, we define the primordial power spectrum as 
(£(ki)£(k2)) = (27r) 3 P^(A;i)(5 3 (ki + IC2) • The average power per logarithmic interval in Fourier 
space is given by V^(k) = j0?z P({k) , and is roughly independent of wavenumber k. The primor- 
dial bispectrum is given by (C(ki)((k 2 )((k3)) = (2ir) 3 B(ki, k 2 ,k 3 )5 3 (k 1 +k 2 + k 3 ) . The ampli- 
tude of the bispectrum relative to the power spectrum is commonly parameterized in terms of the 
non-linearity parameter, /nl, defined such that B(k±, fo, &3)=(6/5)/nl \P{ki)P{k2) + 2 perms] . 

As an example we will consider the primordial curvature perturbation produced by the decay 
of two scalar fields a and b. Without loss of generality, we assume that the curvaton a decays 
first when H = T a followed by the decay of the curvaton b when H = Tb, where r^, < r a . 

During inflation the curvatons acquire an almost scale- invariant spectrum of perturbations on 
super-Hubble scales, Vga* = P&b, = -ff*/(2vr) 2 , where is the Hubble rate at Hubble exit. We 
assume that the field perturbations at Hubble exit have Gaussian and independent distributions, 
consistent with weakly-coupled isocurvature fluctuations of light scalar fields. When the Hubble 
rate drops below the mass of each curvaton field, the field begins to oscillate. The local value 
of each curvaton can evolve between Hubble-exit during inflation and the beginning of the field 
oscillations. We parameterize this evolution by functions g a and g^, but we assume the two fields 
remain decoupled so that their perturbations remain uncorrelated. Thus at the beginning of 
curvaton oscillations, the curvaton fields have values a osc = 9a{o-*) and 6 Q sc = 9b(b*) ■ 

We can define non-linear perturbation ( a and for each curvaton analogous to the total 
perturbation (fTJ. They become constant on scales larger than the Hubble scale once each 
curvaton starts oscillating, and while we can neglect energy transfer to the radiation, r a < H. 
We find that the power spectra of Ca and Cb when the curvatons start to oscillate are related by 

= P 2 P( a , where (3 = ^Jg' w ^h g'^ = dg^/db^ and similar for g' a . In this talk, we assume linear 
evolution between Hubble exit and the beginning of curvaton oscillations, so that (3 reduces to 
the ratio of the background curvaton field values at Hubble exit, (3 = a*/6* . 

3 Full non-linear equations 

We will estimate the primordial density perturbation produced by the decay of two curvaton 
fields some time after inflation has ended using the sudden-decay approximation, generalizing 
the non-linear analysis of Ref.^ to the case of two curvatons. Before the decays we have three 
non-interacting fluids with barotropic equations of state and hence three curvature perturbations 
(PQ) which are constant on large scales 

C 7 = <5iV + im(g) , Ca = ^+|lnfe), C6 = <5iV+iln(|) . (2) 

On the spatial hypersurface where H = T a , there is an abrupt jump in the overall equation 
of state due to the sudden decay of the first curvaton into radiation, but the total energy density 
is continuous 

P71I + Pbi = P 70 i + Pal + Pbi ■ (3) 

Here p 7l 1 is the radiation energy density immediately after the first curvaton decay, pbi is density 
of the second curvaton at the first curvaton decay, and p l0 \ and p a \ are densities of pre-existing 
radiation the first curvaton just before the first decay where a is converted into 71. Note that 
the decay hypersurface is a uniform-density hypersurface and thus, from Eq. ([!]), the perturbed 



expansion on this hypersurface is <5iV = d, where d denotes the total curvature perturbation 
at the first-decay hypersurface. Hence employing Eq. ([2]) we have on this hypersurface 

Pioi =Ay ie 4 ^o-?i), ^ = ^4(0* -CO j p al = p al e 3 ^-^ , p bl = p bl e 3 ^-^ - (4) 

On a uniform-density hypersurface the total energy density is homogeneous. Therefore an 
infinitesimal time before the first decay we have p 7o i + p a \ + Pbi = Pi- Substituting here Eqs. 
and dividing by the total energy density pi we end up with (assuming flat 17 = 1 universe) 

^7oie 4(?7 °- Cl) + Q a ie 3 ^-^ + n bl e 3 ^-^ = 1 , (5) 

where f2 7o i = p 7o i/pi, Q a i = Pal/ Pi, and Q b i = Phil Pi are the energy density parameters of the 
pre-existing radiation, the first curvaton, and the second curvaton at the first curvaton decay, 
respectively. An infinitesimal time after the first curvaton has decayed into radiation we have 

O^e^-CO + n bl e 3 ^-^ = 1 . (6) 

The total curvature perturbation £i is the perturbation on the decay hypersurface and thus is 
continuous between the two phases. However, the energy density of radiation changes abruptly 
at the decay time and the radiation curvature perturbation is discontinuous, C 7o ^ £ 7l . Eqs. ([5]) 
and ([6]) can be used also for the second curvaton decay with the obvious substitutions ft b i — ► 0, 
index a — > b, 1 — > 2, — ► 1. Results for single curvaton follow by setting the density parameters 
of b zero. In Ref. ^ we have checked in single-curvaton case that a fully nonlinear numerical 
solution for a gradual decay agrees well with the analytical solution obtainable from ([5]) and ([6]) . 

If one is interested in perturbative results, one can Taylor expand the exponential functions, 
e x = 1 + x + x 2 j 2 + . . ., and substitute the expansions of the curvature perturbations up to any 
desired order, and then solve order by order for the final curvature perturbation Qi = C72 as a 
function of C 70 , Ca> and Cb- Indeed in the pure curvaton model the pre-existing radiation from 
inflaton decay is homogeneous £ 7o = 0, which we assume from now on. 

4 Results: Primordial power spectrum and /nl 

The first order curvature transfer efficiency parameters^ can be given in terms of three ratios 
which specify the relative densities at the first decay f a i = 3Sl a i / [4O 70 i + 3Q a i + 3fi&i], fti = 
3£l b \/ [4S7 7o i + 3ft a i + 3S7fti] and at the second decay fa = 3Q 2/ [4^ 7l 2 + 3f2{,2]- 

Solving Eqs. ([5]) and ([6]) up to first order, we obtain the primordial power spectrum af- 
ter the second curvaton decay ^ = [r 2 + f3 2 r 2 ] , where (3 2 gives the ratio between the 
initial power in the curvaton a and that in the curvaton b. Here the total first order cur- 
vaton perturbation transfer efficiencies are r a = [(1 — /b2)(3+/ a i)/ai] / [3(1 — hi) + fai] and 
r b = [(l-/&i)/62(3+/oi) + fbifal] I [3(l-/6i) + fed]- Already this first order result is much 
more cumbersome than in the single curvaton caseQ 

In Ref.^1 we solve Eqs. ([5]) and © step by step up to second order for the first and second 
decay, and finally give an expression for the primordial non-linearity parameter /nl after the 
second decay. It is a very complicated function of (3, f a i, f^i and fa- In 3 panels in Fig. Q] we 
show it for (3 = 00, 0, 1 and for three choices of f^i (with different line styles), (a) (5 — > 00: in this 
case the first curvaton a is almost homogeneous. Then the result is very close to the standard 

6 Note that our result differs from that presented recently by Choi and Gong^due to the presence of the extra 
term fti which they (implicitly) assumed zero. The term fbi arises due to the difference between the uniform 
total density hypersurface and the uniform radiation density hypersurface when curvaton a decays, if the density 
of the curvaton b is not negligible. When fibi = we we recover the simpler result® r a = (1 — fb2)fai and r;, = fb2- 

c If the other curvaton, say 6, is absent, then we recover the standard single-curvaton first-order result Pq = 
r a P Ca with r a = f al = 30a/ [40^1 + 3fi„i] = 30a/ [4 - n„i]. 
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Figure 1: /nl(/m, fai) for f) = oo (left panel), (3 = 1 (middle panel), and (3 = (right panel). The energy density 
ratio of the first curvaton at the first decay is set to fbi = (thick solid lines) or fbi = /m/2 (dotted lines) or 
fti = (1 + /ai/3)/t2 (thin solid lines) which corresponds to the simultaneous decay of a and b. Contours of equal 
/nl are shown, from left to right (or from bottom to top), for /nl 



single-curvaton result (in Fig. fa = 0); f^ ee {f) = jj — § — with / = fa. However, if the first 
curvaton density at the first decay is non-negligible, it modifies the result slightly. Large /nl is 
obtained if b is highly subdominant (fa *C 1) at its decay time, (b) (3 = 0: both curvatons carry 
equal amount of perturbation. The result is what one would expect from single-curvaton studies. 
Large /nl is obtained if and only if both a and b are highly subdominant (/ al <C 1, fa <C 1) at 
their decay time, (c) (3 = 0: the second curvaton is homogeneous. In this case we can obtain 
a large /nl even if both the curvatons are dominant (/ al ~ 1, fa ~ 1, upper right corner in 
Fig.) at their decay time. The inhomogeneous radiation produced by the first curvaton decay 
is diluted by the decay of the second homogeneous curvaton. We expect this result to hold also 
for more than two curvatons if the last decaying curvaton is almost homogeneous. 

In all cases we find /nl > — 5/4, which seems to be a robust lower bound in single and 
multi-field curvaton models in which the curvaton field perturbations are themselves Gaussian. 
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